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(N 

Using a general relativistic exact model for spherical structures in a cosmological background, we 
^ I have calculated the test particle geodesies within the structure for different masses in order to obtain 

Qthe velocity profile of stars or galaxies. Defining a Newtonian mass based on the classical dynamical 
' relations, it turns out that the Misner-Sharp quasi-local mass is almost equal to the Newtonian one. 

' This, however, is not the case for other general relativistic quasi-local mass definitions, which can 

be much smaller than the mass definition based on the classical dynamics. Therefore, based on the 
rotation curve, we are not in a position to relate a unique mass to a cosmological structure within 
general relativity even in cases of very weak gravity. 



X 



PACS numbers; 



I. INTRODUCTION 

On ' How good is gravity described by Newtonian gravitation versus Einstein's general relativity (GR) at scales larger 
C than solar system? Within the solar system we know already what the effects of general relativity are. It is also 

widely accepted that at larger scales and in regions far from strong gravity fields, such as within galaxies or clusters 
C\| , of galaxies and even at the cosmological scales as used in different cosmological simulations, the Newtonian gravity 
' is valid to a large extend. There are, however, contrary claims that because of the non-linear character of Einstein 
equations this may not be true. Wald et.al. have faced this problem for cosmological cases in an elegant albeit general 
' , [ • approximative way There are still some facts related to the non-linearity of the equations which cry for exact 
J> ■ models within GR to be compared with the corresponding Newtonian ones. Any final and definite conclusion about 
the validity of the Newtonian approximation depends on the results of such exact models and the comparison to their 
corresponding Newtonian approximations. As a contrary example of the validity of Newtonian approximation at large 
, scales, take for exarnple the mass definition in general relativity which is not unique and may differ substantially from 
the Newtonian one [lOl|. Another example is the difference between the unique isolated horizon of the Schwarzschild 
metric and the corresponding multiple dynamical horizons in a cosmological black hole where a "negligible" cosmic 
mass density of the^order of 10~^^gr/cm'^ outside the horizon makes this global difference to the vacuum solution of 
the Schwarzschild a difference not explained within the Newtonian approximation. These phenomena could only 
be explained by taking into account the effect of the non-linearity of GR in the presence of an extra week gravity 
field and by having exact solutions to be compared with the corresponding approximations. 

We intend to construct a simple exact astrophysical GR model of a structure within a FRW universe, study its 
geodesies and the corresponding mass definitions, and compare the result with the Newtonian counterparts to see 
if and to what extent the familiar claims of the applicability of the Newtonian gravity is correct. Our concern 
is neither dark matter nor dark energy but to what extent and in which sense the linearized Einstein equations are valid. 
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There have been different endeavors to construct solutions of Einstein equations representing a coUapsing central 
mass. We are interested in an exact solution of Einstein equations representing a cosmological structure as an overden- 
sity region within a cosmological background which we assume to be asymptotically FRW. There are not much viable 
exact models representing such a structure. We will use a recent analytical model proposed by the authors based on an 
inhomogeneous cosmological LTB model to construct such an asymptotically FRW universe including an overdensity 
region evolving to a black hole different from that of the Schwarzschild [3] . Our aim is to understand the notion of mass 
related to different definitions of quasi-local masses for such a cosmological structure within a FRW universe, based 
on the behavior of geodesies and the corresponding rotation curve, and compare them to the classical Newtonian mass. 

The Arnowitt-Dcscr-Misner mass [3] (MjijjM) at the spatial infinity and the Bondi-Sachs mass 4] (Mbs) at 
the null infinity of asymtotically Minkowskian space-times are among the early definitions of mass which are not 
extendable to non-static and non-asymtotically fiat space-times, where only a so-called quasi-local definition is possible. 
Hawking defined a quasi local mass which has various desirable properties: it is zero for a metric sphere in a fiat 
spacetime, gives the correct mass for the Schwarzschild solution on a metric sphere, and tends to the Bondi-Sachs 
mass asymptotically in a static asymptotically flat spacetime. However, it is non-zero for generic 2-surfaces in a 
fiat spacetime. Another attempt was due to Misner and Sharp 5] (Mms) with a well- understood Newtonian limit. 
One of the most promising quasi- local definitions, however, seems to be the one proposed by Brown- York (Mby) 
0. Motivated by the Hamiltonian formulation of general relativity, they found an interesting local quantity from 
which the definition of quasi-local mass was extracted. Their definition depends, however, on the choice of the gauge 
along the 3-dimensional spacelike slices. It has the right asymptotic behavior but is not positive in general. Using 
some geometric considerations, two new quasi-local mass definitions were put forward by Liu-Yau {M^y) and Epp 
(MeppM) which are gauge independent, and always positive. These quasi-local masses differ for a cosmological black 
hole as discussed in . We will see that some of these mass definitions are equal to each other in our spherically 
symmetric case. 

There has recently been some attempts to study general relativistic rotation curve (see ,11] and pJl). The model 
discussed in these papers, however, is just based on the simple Schwarzschild metric, dust collapse, or a simple axially 
symmetric space-time with a singularity at the central plane, without looking at a model structure within a FRW 
universe. To achieve this, we start in section II with a brief introduction to the LTB metric, followed by a definition 
of quasi-local masses in section HI. In section IV we introduce then our model for a test particle motion around a 
dynamical cosmological black hole within a FRW universe. The rotation curve for a more realistic model is obtained 
in section V. We then conclude in section VI. Throughout the paper we assume SirG — c — 1. 



II. LTB MODEL OF A STRUCTURE 



The LTB metric in synchronous coordinates is written as 



c/s2 ^ ?-—dr^ - i?(r, tfdn'^, (1) 

1 + .f{r) 

representing a pressure-less perfect fiuid satisfying 

Here dot and prime denote partial derivatives with respect to the parameters t and r respectively. The angular 
distance R, depending on the value of /, is then given by 

M 

R = — ^(1 - cosT]{r,t)), 

(_f)3/2 

V-smv= ^ ■'^^ {t-tb{r)), (3) 



for / < 0, and 



R^ilM)Ht-h)K (4) 
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for / = 0, and 

M 

R = —{cQshri{r,t) - 1), 

J3/2 

sinh?7 - 7? = -^(t - tf,(r)), (5) 

for / > 0. 

The metric is covariant under the rescahng r — > f(r). Therefore, one can fix one of the three free functions of the 
metric, i.e. /(r), or M{r). The function M{r) corresponds to the Misner-Sharp mass in general relativity, as 

shown in the general case of spherically symmetric solutions of Einstein equations [H. The r dependence of the bang 
time tb{r) corresponds to a non-simultaneous big-bang or big-crunch singularity. 

There are two generic singularities of this metric, where the Kretschmann and Ricci scalars become infinite: the 
shell focusing singularity at R{t,r) = 0, and the shell crossing one at R'{t,r) — 0. However, there may occur that 
in the case of R{t, r) = the density p = -^2^ and the term ^ remain finite. In this case the Kretschmann scalar 
remains finite and there is no shell focusing singularity. Similarly, if in the case of vanishing R' the term is finite, 
then the density remains finite and there is no shell crossing singularity either (see for more detail). 



III. QUASI-LOCAL MASSES 

There has been many attempts to define a quasi-local mass (QLM)or energy in a dynamical setting, where no 
global definition is possible. It has been turned out that there is no unique definition of a QLM. In fact more than 
10 definitions have already been proposed in the literature. The difference between some of the proposed quasi-local 
mass definitions has been studied in jici| . We have seen in that in the case of spherically symmetric structures the 
Hawking mass is equal to the Misner-sharp once. Here we give a brief introduction to the most familiar quasi-local 
mass concepts. The Brown- York mass for the 2-boundary specified by r = costant and t — constant (MBY.r) in an 
asymptotically FRW solution is given by 



Mby.v = -i? VTT7 + (RVT+f) \frw ■ (6) 

If we specify the 2-boundary by R{r, t) — costant and t = constant with R being the physical radius, then the 
Brown- York mass {Mby.r)) is given by 



2M / 2M 

Mby.r = -R^ ^ + IT ^ V ^ ^ IT^ '^'^'^ ' 

The Liu-Ya and Epp masses are equal in our spherically symmetric case and is given by 



2M / 2M 

Mly = Mepp = -R\ll + — + 1 + — ) \frw ■ (8) 



We therefore concentrate on three mass definitions: the Misner-Sharp one which is equal to that of Hawking; The 
Epp mass being equal to Liu-Yau and Brown- York with constant physical radius (i? = constant); and finally the 
Brown- York for constant comoving radius (r = constant) . 

IV. CONSTRUCTING THE MODEL 



The toy model we are going to construct and study should describe a simple model of a spherically symmetric 
mass condensation within a FRW universe as an exact solution of Einstein equations with a pressure-less ideal fiuid. 
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Therefore we will choose a density profile reflecting an overdensity at the center and almost constant density far 
from the center as expected for a FRW universe. Within such a model structure we then study timelike geodesies to 
extract information about the rotation curve in such a dynamical setting. This is achieved by specifying the three LTB 
functions tb{r), f{r), and M{r). Assuming an expanding universe, this cosmic LTB model structure starts expanding 
with the universe before its expansion decouple from the universe and a collapsing phase starts. There are different 
ways to specify the LTB functions depending to our needs or our methodology. Once the solution is fixed we may 
study the timelike geodesies to extract the rotation curve within the structure. Now, the geodesic equations for the 
LTB metric are given by 



and 



d'r _ m+J) , d<l> , d0 , i^'^"-Wwl dr , R' dt dr 
d>? - — R^NnW (^) +id>)\ ^2 (^) - 2;r7(^)(^). (9) 

d^6 . ,^^,d(l),-2 „R'dr.,d9. ^R'dt.,d9. 

d>? = '''^'M^ - 'r^Tx^^Tx^ - 'r^Tx^^Tx^^ ^''^ 

fjt = r, cos{0) , d9 d(j) _ 2:^/^w^\ 2:(^C^^C^^ nn 
dX^ sm{e)^ dX'^dx' R^dx'^dX' R^dX'^dx' ^ ' 

Using these geodesic equations, we arc able to find dynamical properties of a test particle within this structure. 
In addition, we may use any general relativistic mass definition to calculate the mass of the structure. It then 
remains to find the corresponding Newtonian mass to be compared with the general relativistic ones. Hence, we use 
the Newtonian relation to adhere a mass to the structure leading to the circular velocity reflected in the geodesic 
equations. The Newtonian mass is obtained according to 

= (13) 

Given that our cosmic structure is dynamic we choose the initial conditions such that the structure is in its late 
expanding phase. Therefore, we may expect to have quasi-circular geodesies. These arc defined either as those having 
a vanishing radial velocity and acceleration at the initial conditions, or have an almost constant radius within the 
numerical prcicision for a finite angular displacement. We have checked both procedure leading to the same numerical 
result. Note that we have assumed R to be the physical radius with A its afflne parameter and the acceleration given 
by 

dt- (^)3 ^ ^ 

The terms ^) 3X' 31 ^^'^ calculated from the geodesic equations. 



V. MASS AND ROTATION CURVE FOR THE GENERAL RELATIVISTIC COSMIC STRUCTURE 

We report the result in three stages. First we apply our method to the Schwarzschild case just to test our code. 
We then go over to a toy model structure and finally apply the method to a more realistic case with a NFW density 
profile. 



A. rotation curve in the vacuum case of SchwEirzschild 



Noting that the Schwarzschild space-time is a particular case of LTB space-time, we choose the three LTB functions 
in the following way: 



fir) = 0, 



(15) 
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M(r) = m, 
tb{r) = r. 



(16) 
(17) 



The resulting Schwarzschild metric is then given by 



(|m(r-i))^ 



dr^ + (^m(r - tf)i{de^ + sin{efd(j)^). 



(18) 



This is the metric of the Schwarzschild space-time in a synchronous frame. By studying circular geodesies having 
vanishing radial acceleration and velocity, we obtain the rotation curve as given in the Fig.(IT]). Here the mass is the 
unique Schwarzschild one. Obviously the resulting rotation curve in the region 2gL << 1 is the same as what we 
expect from the classical mechanics, showing the validity of our numerical method and at the same time the equality 
of the Newtonian mass and the the general relativistic Misner-Sharp mass. 



M= 1.069 10'"M 



- Newtonian curve 



C 
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FIG. 1: Rotation curve for Schwarzschild metric 



A Toy Model 



The toy model is constructed such that it represent a cosmic structure within a FRW universe as an exact solution 
of the Einstein equations. We fix the three LTB functions in the following wayQ: 



M{r) 



„3/2 



tb{r) = 0, 



(19) 



The model is a typical cosmic structure showing the formation of a central black hole with a distinct even horizon 
and apparent horizon, and a void before entering the asymptotic FRW region. Fig.© shows the density profile of the 
toy model and its corresponding rotation curve. 

We now associate classical and relativistic masses to this curve. Whereas the classical Newtonian mass is unique, 
in general relativity we are free to choose any of the many definitions put forward in the literature (see for example 
the discussion in [ifll)- We concentrate on those quasi-local mass definitions put forward in section III. The result is 
shown in Fig. ([3]), where the ratio of general relativistic masses to the classical Newtonian one is also depicted. As we 
can see from those figures, the ratio of the Mms rnass over the classical mass is almost 1. In other relativistic cases, 
however, the ratio may even decrease down to about 0.2. Therefore, given a rotation curve, the interpretation of it 
based on general relativistic dynamics may lead to a corresponding mass of the structure which maybe much less 
than the Newtonian one depending on the definition we choose. 



6 



FRWDensity \{y 



5 10 50 100 500 

Physical Radius (kpc) 



Physical Radius ( kpc ) 



FIG. 2: Density profile for a toy model of galaxy and rotation curve related to it 
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FIG. 3: The Quasi-local masses for a toy model of galaxy 



C. NFW model 



Now, we try a model structure with a NFW density profile [l5|. We use the this density profile to find the physical 
radius through 



R{U,m) = / 
Jo 



37rp{t, x) 



-dx. 



(20) 



Assuming then the density profile at two different times ^1,^2, i-e. p{ti,r) > p{t2,r), we obtain the corresponding 
physical radii R{t2,r) < R {ti,r), leading to 



p^{r) = Pcr;t{U){{5cMBe ^ " 6i)e + 1), 



(21) 



P/W = Pc«t(t/)(( 



5c 



"^«/^' + 1). 



(22) 



Choosing the third function in the LTB model, the two physical radii results in the numerical solution of our the LTB 
model structure. Fig. ([4]) shows the resulting density profile and the related rotation curve of the NFW model for a 
typical galaxy mass. 



We are now in a position to calculate different masses. The results are shown in Fig.®. As in the case of the 
toy model, we realize that the Misner-Sharp mass is almost identical to the Newtonian mass. There is however 
substantial difference to the general relativistic masses. The ratio of these quasi-local masses to the Newtonian one 
again may be as small as 0,2. 
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FIG. 4: Density profile for the NFW model of a galaxy 
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FIG. 5: Ratio of GR masses and classical masses for the NFW model of a galaxy 

VI. CONCLUSION 

Our study of exact general relativistic cosmic structures and their effect on the motion of test particles, specially 
the rotation curve, shows a variety of interesting departures from the Newtonian dynamics. The classical Newtonian 
mass related to a sample rotation curve is calculated and compared to different corresponding general relativistic 
mass definitions. While the Misner-Sharp and the Hawking mass are almost identical to the Newtonian one, the 
other masses such as Brown- York, Liu-Yau, and Epp are smaller in a wide range of distances from the center of 
the structure. These general relativistic masses corresponding to the sample rotation curve may even be less than 
0.2 of the corresponding Newtonian mass. We conclude that adhering a mass to a structure with a definite density 
profile based on the rotation curve is conceptually controversial within general relativity for astrophysical mass scales. 
Whether this has any effect to the problem of dark matter has yet to be discussed. We have, however, to conclude 
that there is no unique mass definition within general relativity based on the rotation curve and for large scales. If 
there is a need to define a mass as an integrated quantity at all in such astrophysical and cosmological scales is then 
a matter of dispute. 
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